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Abstract 

Monte Carlo simulations have been performed to determine the excess energy and the equation 
of state of fee solids with Sutherland potentials for wide ranges of temperatures, densities and 
effective potential ranges. The same quantities have been determined within a perturbative scheme 
by means of two procedures: i) Monte Carlo simulations performed on the reference hard-sphere 
system and ii) second order Barker-Henderson perturbation theory. The aim was twofold: on the 
one hand, to test the capability of the 'exact' MC-perturbation theory of reproducing the direct 
MC simulations and, on the other hand, the reliability of the Barker-Henderson perturbation 
theory, as compared with direct MC simulations and MC-perturbation theory, to determine the 
thermodynamic properties of these solids depending on temperature, density and potential range. 
We have found that the simulation data for the excess energy obtained from the two procedures 
are in close agreement with each other. For the equation of state, the results from the MC- 
perturbation procedure also agree well with direct MC simulations except for very low temperatures 
and extremely short-ranged potentials. Regarding the Barker- Henderson perturbation theory, we 
have found that, surprisingly, the first-order approximation is in closer agreement with simulations 
than the second-order one. 

PACS numbers: 05.10.Ln, 05.70.Ce, 64.10. +h 
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I. INTRODUCTION 



The Sutherland potential is of the form 



u(r) = 




) 



7 



r < a 



r > o 



(1) 



were parameters a and — e account for diameter of the particles and the maximum potential 
depth, respectively, and the exponent 7 determines the effective range of the potential, so 
that the higher is the value of 7 the shorter is the effective range. 

The shape of the potential function (JTJ) resembles the shape of other widely used poten- 
tials, such as the Lennard- Jones potential, for 7 = 6, the Girifalco potential for Cqq, for 
7 « 12, 1 and the Mie potential, widely used to describe pure fluids and solutions as well as 
solids, for a wide range of values of 7 depending on the real system to be modelled.- 2 The 
advantage of a potential of the form ([1]) with respect to these others is that in the former 
the diameter of the particles is well defined. The simple mathematical form of the potential 
(CQ), together with its flexibility to reproduce a wide variety of intermolecular interactions, 
has led recently to a number of applications for modelling interactions in complex fluids 3,4 
and nanocomoposite materials. 5 

In this situation, one might reasonably expect this kind of systems would have been 
extensively studied bot from theory and from computer simulation. However, computer 
simulations on this system have been scarcely reported. For the fluid, quite extensive simu- 
lation data for the excess energy and the compressibility factor have been reported in ref.~ 
and more recently in^ for several values of 7 and a wide range of densities and temperatures, 
particularly in the latter reference. Simulation data on the radial distribution function for 
the high density fluid may be found in ref.- for several temperatures and effective ranges. 
Simulation data for liquid-vapour coexistence densities and critical parameters for several 
values of 7 were reported in ref.- and triple point parameters in ref.—. 

Regarding theory, the amount of research devoted to this kind of fluids is equally 
scarce^ 1 ^ 1 ^ 3 -, although recently we reported^ a quite complete analysis of the performance 
of the the Barker-Henderson perturbation theoryi^ 1 ^ for the thermodynamic properties of 
this kind of fluids, and still more recently^ we performed a similar analysis using the first- 
order mean spherical approximation.— »22 
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In contrast with fluids, crystalline solids interacting by means of model potentials have 
been much less frequently studied, with the possible exception of the hard-sphere solid, and 
to the best of our knowledge the behaviour of solids with Sutherland potentials has not 
been studied before from theory nor computer simulation. This is the aim of this paper, 
for which we have performed extensive computer simulations for the thermodynamic prop- 
erties of Sutherland fee solids for several values of 7, temperatures and densities. These are 
described in the next Section. Section III, summarizes the foundations of two versions of 
the perturbation theory 1) an 'exact' first-order perturbation theory based on simulations 
performed on the reference hard-sphere fee solid and the Barker-Henderson perturbation 
theory. These two theories are compared in Section IV with each other and with the sim- 
ulation data reported in Section II, and the concluding remarks are presented in the same 
Section. 

II. MONTE CARLO SIMULATIONS FOR FCC SOLIDS WITH SUTHERLAND 
POTENTIALS 

We have performed NVT Monte Carlo simulations for fee solids with potentials of 
the form (CD) with 7 = 6, 12, 18, and 36, reduced temperatures T* = 0.6, 0.8, 1.0, 1.5, 2.0 and 
3.0, and reduced densities p* in the range 0.90 — 1.30 with step 0.05 this covering most of 
the density range for the solid and some of the high density fluid. The potential cut-off 
distance was fixed at r c = 3a. Most of the systems consisted in N = 500 particles but 
some simulations were carried out with 2048 particles to test the influence of the size of the 
system. The particles were placed in a cubic box with periodic boundary conditions in a fee 
configuration. The system was equilibrated for 2 x 10 4 cycles, each cycle consisting of an 
attempted move per particle and then the thermodynamic and structural properties were 
measured over the next 5 x 10 4 cycles, with partial averages every 500 cycles from which the 
statistical errors were determined as the standard deviation. Acceptance ratio was fixed at 
around 50%. The excess energy U E was determined from the energy equation in the form 




(2) 



The compressibility factor Z was determined from the expression 
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whichh results from the combination of eq. ([2]) with the virial theorem for the Sutherland 
potential. 

Results are listed in Tables Hl lIVl in which corrections to account for the effect of the 
truncation of the potential have not been included. For the Sutherland potential with cut- 
off distance r c = 3<r, the corrections for the compressibility factor and the excess energy 
are 

AZ = -tt T - P- (4) 

3 3-7T* K ' 

and 

AU E 3 3 ~^ 

respectively. 

III. PERTURBATION THEORY 

In perturbation theories, the intermolecular potential is split in the form 

u(r) = «o(r) +wi(r), (6) 

where Uq(t) and U\(r) accounts essentially for the short-range and long-range contributions, 
respectively. The former is mainly due to the repulsive forces and the latter to the attractive 
ones. For the Sutherland potential (pQ), the obvious choice for these contributions is 



, r < a 

u (r)={ ~ , (7) 

r > a 

and 



r < a 

«i(r)H /x7 , (8) 

respectively. Note that ([7]) is the hard sphere (HS) potential. 
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At high densities the thermodynamic and structural properties of a system are deter- 
mined mainly by the first of these contributions, whence arise perturbation theories. These 
theories consider that the properties of the system at high densities are given by those of 
a reference system, one consisting of particles interacting by means of the potential Uq(t) 
and providing the main contribution, plus a minor contribution due to M 1 (r) and considered 
as a perturbation of the former. As most of the thermodynamic systems with spherically- 
symmetric potentials at high temperatures behave much like a hard-sphere system, the latter 
is the obvious choice for the reference system. Therefore, it is natural to express the ther- 
modynamic properties as a series expansion in terms of 1/T*, the inverse of the reduced 
temperature, with the zero-order term corresponding to the contribution of the reference 
system and the remaining terms accounting for the contribution of the perturbation. In this 
situation, the Helmholtz free energy, the radial distribution function, the excess energy, and 
the compressibility factor can be expressed in the form 



oo 

V— - - (9) 



Nk B T ^ Nk B T T 



n=0 



9{r) 



x ^ 

J>n(rW, (10) 



n=0 

tjE °° TJ i 

^L = Y—— (11) 

n=0 

oo .. 

Z = J2 Z ^ ( 12 ) 

n=0 

where subscripts '0' correspond to the contributions due to the reference system which, for 
systems with the Sutherland potential, is the hard-sphere system. 



A. 'Exact' First-order Perturbation Theory 

Before analysing the performance of the Barker-Henderson perturbation theory for 
Sutherland solids, it is worth analysing the capability of an 'exact' perturbation theory 
for obtaining the thermodynamic properties of these solids. By 'exact' by mean a pertur- 
bation theory with the terms in expansions (|T0|) - (fl2|) obtained from computer simulations 
performed on the hard-sphere reference system, thus avoiding any theoretical approximation. 
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The procedure to determine the zero- and first-order terms in the expansion ({TO]) by 
means of computer simulation was developed by Smith et al.— quite a long time ago. They 
found 

0b(r< + Ar/2) = - (13) 
27iNp(rf +1 - rf) 

and 

^i(r< + Ar/2) = * (14) 

2ti Np (rf +1 - rf) 

In the preceding expressions Ni is the number of intermolecular distances in the range 
( r i, r i+i), with Ar = r i+1 — rj << a i = 0,1, . . ., angular brackets mean an average, subscript 
mean that the averages are performed in the reference HS system, and u*(r) = u\{r)/e. 

Introducing the r.d.f. (TI01) . truncated at first order, into the energy equation ([2]) provides 
the first- and second-order contributions to the excess energy in the form 

U x = 2N<K P Y,9o{n)r*u 1 (r i )Ar (15) 

i 

and 

U 2 = 2N?rpJ29i(ri)r*u 1 (r i )Ar (16) 

i 

respectively. 

Introducing in turn the excess energy (llip . truncated at first order, into the virial equation 
of sate ([3]), provides the zero- and first-order contributions to the compressibility factor. 
They are 

Z = l + ^p^g {a) (17) 

and 

z, = l^sM + (18) 

respectively. 
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In principle, the procedure could be generalised to obtain higher-order terms, but the 
computational effort increases with the order of the perturbative term and becomes imprac- 
tical. 

This procedure, that has been denoted Monte Carlo pertubation theory (MC-P), 21 was 
applied to obtain the first terms in the perturbative expansions (ITU|) - flT2"|) for square- well 
fluids, first by Smith et al.— and more recently by two of us.— i22j2I Very recently we have 
used successfully the same procedure to obtain the thermodynamic and structural properties 
of fluids with Sutherland potential.-. To the best of our knowledge, the procedures has never 
been applied to crystalline solids. 

The explicit expressions of eqs. (}TD]) - (fT2|) . truncated at the above-mentioned level, are as 
follows 

9 0) = 9o 0) + 9i (x) (19) 
where x = r/a is the reduced distance, 

^ = El + Ei± (20) 

iVe iVe NeT* V ' 

and 

Z = Z + ±Z X (21) 

We have performed Monte Carlo simulations in the NVT ensemble of the reference HS 
fluid to obtain the perturbative terms in expressions (1191) - 02 II) . To this end, we considered a 
system consisting in 500 hard spheres in an placed in a fee configuration within a cubic box 
with periodic boundary conditions. At each density, the system was allowed to equilibrate 
for 5 x 10 4 cycles and the averages involved in expressions (fl3|) and (fl4|) were calculated from 
the next 10 6 cycles. We chose Ar = 0.005 and the acceptance ratio was settled to about 
50%. From these data, we have calculated the values of go(x) an d giix) from eqs. (TT3l) and 
dEJ), respectively, and subsequently Ui, U 2 , Z , and Zi from eqs. f|T5l) - ffTSl) . 
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B. Barker-Henderson perturbation theory 



The Barker-Henderson perturbation theory (BH)^^ provides expressions for the flrst- 
and second-order terms in the expansion ([9]). The first-order term is given by 



Fx 



POO 

2np / u\(r)go(r)r 2 dr (22) 
Jo 



Nk B T 

The second-order term in the macroscopic compressibility approximation (mc) is 



Nk B T 

and in the local compressibility approximation (lc) 



-<K P k B T j™ [ul(r)) 2 [^j g (r)r 2 dr (23) 



F2 --*pk B T ri<(r)?( dlP9 : {r)y ) r 2 dr (24) 



o 



Nk B T r Jo 1 \ dp 

The corresponding contributions to the excess energy can be readily obtained from the 

thermodynamic relationship 



(25) 
J v 



(26) 



= 2 \ d{F/Nk B T) 
Ne [ dT* 

and those for the compressibility factor from 

z = P" V _ \ d(F n /Nk B T) 
n Nk B T P [ dp 

Using suitable expressions for the r.d.f. go(r) and the equation of state Z of the HS 
reference fluid, Henderson and coworkers^ 1 ^ 1 ^ 1 ^ 1 ^'^, and many others, applied successfully 
the BH theory to a variety of model fluids. 

The Barker-Henderson theory has been sometimes applied to crystalline solids.— i2&2iL2£ 
To apply this theory, accurate expressions for the compressibility factor Zq and the r.d.f. 
go(r') of the hard-sphere reference fee solid are needed. The equation of state of a fee hard- 
sphere solid is well reproduced by the Hall 31 equation 

Z = 3/a + 2.557696 + 0.1253077/3 + 0.1762393/5 2 - 1.053308/5 3 + 2.818621/3 4 

-2.921934/3 5 + 1.118413/3 6 (27) 

where (3 = 4a/ (1 + a) and a = p / p — 1, in which p is the close-packing density. 



Regarding the radial distribution function of the fee hard-sphere solid, accurate 
parametrizations are available.— MMi2k One of the most frequently used is 3 - 3 - 

<7o (x) = — exp [— (x — x\) — W 2 4 (x — Xi) ] 

where 77 = (7r/6)p* is the packing fraction and X\ is the position of the first peak of the 
r.d.f., which was determined 3 ^* 3 " 3 - from the simulation data for several packing fractions. 
Functions Wi, W2, and W are given by 36 



For 77 < 0.55: 



/3 

W 1 = V3W 2 = 5 V Q552 exp [10.49375 (77 - 0.52)] (29) 



[0.23601 (77 - r/ ) 2 - 21.4395 (77 - r] f] 



w = ~ ; : (30) 



For 0.55 < r] < 0.73: 



1 5522782 

W x = - 2.0302556 exp [5.8331273 (770 - 77)] + 74.873192 (770 - r?) 2 (31) 



_ 0.9559565 - 5.855022 (r] - 77) + 39.74663 (r/ - r/) 2 - 109.62638 (r/ - r]) 3 

W2 — (32) 

?7o - V 



w = ,™ ^ (33) 



1 - 10.589574 (770 - 'if'' 
[0.694 (77o-r/)] ] 

In the preceding expressions 770 is the packing fraction corresponding to the close packing 
density p f° r the fee solid. 

Finally, parameter A in expression (|28p from the condition that the virial theorem 

Z = 1 + 4t7 5o (1) (34) 

must be satisfied. 
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IV. RESULTS AND DISCUSSION 



In Figure [TJ the results for the excess energy from the BH and MC perturbation theories 
are compared with each other and with the MC simulation data listed in Tables HllIV[ and 
a similar comparison is performed in Figure [2] for the compressibility factor. The upper 
limits in the integrals involved in the BH perturbation theory were set at x — 3.0, thus 
making possible a direct comparison with simulation data without need of including in the 
latter the corrections due to the truncation of the potential. Let us first examine the MC 
simulation data. At densities p* < 0.95 all the systems studied underwent a melting tran- 
sition. Moreover, at the lowest temperatures considered the transition showed a tendency, 
clearly seen in the figures, to displace towards higher densities (p* ~ 1.0), and the effect is 
more pronounced the higher is the value of 7. This is a consequence of the clustering of the 
particles of these systems at low temperatures. 

Regarding the MC-P results for the excess energy U E , one can see in Figure [1] that this 
'exact' first-order perturbation theory reproduces very accurately the simulation data for all 
the values of 7, at all temperatures and densities except near melting at low temperatures. 
The same is true for the compressibility factor Z, as shown in Figure [21 except for 7 > 12 
and T* = 0.6, the lowest temperature considered, and perhaps sometimes for T* = 0.8 too. 
The reason for the better performance of the MC-P results for the excess energy than for 
the compressibility factor is, at least partially, because this 'exact' perturbation theory is 
second order in the former quantity whereas it is only first order in the radial distribution 
function and the equation of state. 

Now, let us analyze the performance of the BH perturbation theory. For the excess energy, 
at temperatures T* > 1.5, the contribution of the second-order term is nearly negligible, 
so that the BH theory truncated at first order provides satisfactory results in all cases for 
T* > 1.5 and even for lower temperatures in the case 7 = 6. For 7 = 12 and T < 1.0, second- 
order approximation is needed, with the local compressibility approximation providing grater 
accuracy than the macroscopic compressibility approximation. For 7 > 18 and T < 1.0, 
the second-order approximation fails to provide completely satisfactory results for the low- 
density solid. It is to be remarked that even in the cases where the BH theory departs from 
simulation data, the 'exact' first-order MC-P continues doing a good job. 

As far as the equation of state is concerned, Figure [2] shows that the BH perturbation 
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theory is in complete agreement with the MC-P theory up to densities close to the melting 
density, except perhaps for 7 < 18 and T* = 0.6, with little or no difference between 
first- and second-order approximations, although the local compressibility approximation 
seems to work slightly better at low temperatures than the first-order or the macroscopic 
compressibility approximations. 

From the precedent analysis it seems clear that if we were able of improving the theo- 
retical predictions for the excess energy until reaching the accuracy of the MC-P results, 
and using the Helmholtz free energy route to obtain the equation of state instead of the 
pressure equation, the results obtained from the perturbation theory would probably be in 
complete agreement with simulation data for solids with spherically symmetrical potentials 
even for extremely short-ranged potentials and low temperatures. While this improvement is 
achieved, the Barker-Henderson perturbation theory constitutes an excellent choice in most 
situations. 
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List of Tables 



TABLE I: Simulation data for 7 = 6. Corrections due to the truncation of the potential have not 
been included. The number between parenthesis indicates the error in the last significant digit of 
the compressibility factor. For the excess energy, the statistical error in most cases is beyond the 
third decimal place and when this is the case it is not indicated. 

p* 0.90 0.95 1.00 1.05 1.10 1.15 1.20 1.25 1.30 
T* =0.6 

Z 2.70(2) 3.75(3) 0.53(2) 0.28(2) 0.70(2) 1.66(2) 3.84(2) 8.05(3) 17.41(4) 

U E /Ne -3.382 -3.651 -3.852 -4.149 -4.479 -4.832 -5.208 -5.610 -6.036 

T* =0.8 



Z 4.59(2) 5.96(3) 2.84(2) 3.04(2) 3.75(2) 5.22(2) 7.71(2) 12.42(3) 22.24(4) 

U E /Ne -3.351 -3.630 -3.835 -4.141 -4.473 -4.829 -5.207 -5.609 -6.036 

T* =1.0 



Z 5.85(2) 7.22(2) 4.32(2) 4.72(2) 5.65(2) 7.29(2) 10.05(2) 15.01(3) 25.04(4) 

U E /Ne -3.336 -3.616 -3.826 -4.136 -4.470 -4.827 -5.206 -5.609 -6.036 

T* =1.5 



Z 7.481(2) 9.04(2) 6.23(2) 6.95(2) 8.16(2) 10.07(2) 13.15(2) 18.49(2) 28.91(4) 

U E /Ne -3.319 -3.602 -3.815 -4.129 -4.466 -4.825 -5.205 -5.609 -6.036 

T* =2.0 



Z 8.277(2) 9.77(6) 7.24(2) 8.06(2) 9.41(2) 11.50(2) 14.73(2) 20.21(2) 30.79(4) 

U E /Ne -3.310 -3.593(1) -3.810 -4.126 -4.464 -4.824 -5.204 -5.609 -6.036 

T* =3.0 



Z 9.12(2) 8.18(4) 8.23(1) 9.21(2) 10.74(2) 12.92(2) 16.29(2) 21.92(2) 32.71(4) 

U E /Ne -3.303 -3.529(1) -3.805 -4.123 -4.463 -4.823 -5.204 -5.608 -6.036 
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TABLE II: As in Table Q] for 7 = 12. 



p* 0.90 0.95 1.00 1.05 1.10 1.15 1.20 1.25 1.30 
T* =0.6 

Z 4.54(3) 5.61(3) 5.98(9) 2.52(2) 2.03(3) 1.86(3) 2.36(3) 4.43(3) 10.884(4 

U E /Ne -2.018(1) -2.221(1) -2.432(2) -2.510 -2.760 -3.068 -3.442 -3.892 -4.429 

T* =0.8 



Z 6.06(2) 7.44(3) 5.03(3) 4.61(2) 4.68(2) 5.21(2) 6.54(3) 9.64(3) 17.28(4) 

U E /Ne -1.951(1) -2.167(1) -2.247(1) -2.458 -2.724 -3.046 -3.430 -3.886 -4.428 

T* =1.0 



Z 7.04(2) 8.36(3) 5.87(2) 5.97(2) 6.40(2) 7.31(2) 9.13(3) 12.77(3) 21.195(4) 

U E /Ne -1.919(1) -2.130(1) -2.205(1) -2.431 -2.707 -3.035 -3.423 -3.883 -4.427 

T* =1.5 



Z 8.26(2) 9.88(2) 7.23(2) 7.66(2) 8.61(2) 10.09(2) 12.49(2) 16.96(3) 26.23(4) 

U E /Ne -1.873(1) -2.089 -2.159 -2.396 -2.683 -3.019 -3.414 -3.880 -4.426 

T* =2.0 



Z 8.86(2) 10.59(2) 7.96(2) 8.60(2) 9.72(2) 11.49(2) 14.17(2) 19.10(3) 28.76(4) 

U E /Ne -1.851 -2.070 -2.139 -2.381 -2.671 -3.011 -3.409 -3.878 -4.426 

T* =3.0 



Z 9.46(2) 10.5(1) 8.71(2) 9.53(2) 10.95(2) 12.89(2) 15.93(2) 21.18(3) 31.36(4) 

U E /Ne -1.829 -2.014 -2.119 -2.366 -2.661 -3.004 -3.406 -3.876 -4.425 
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TABLE III: As in Tabled for 7 = 18. 



T* =0.6 

Z 5.411(3) 6.40(3) 6.83(7) 4.17(3) 3.48(4) 3.02(3) 2.68(3) 2.89(4) 6.48(5) 

U E /Ne -1.604(1) -1.778(2) -1.945(3) -1.995(2) -2.191(1) -2.454(1) -2.790 -3.231 -3.822 

T* =0.8 



Z 6.73(2) 8.04(3) 6.81(3) 5.80(2) 5.77(3) 5.95(3) 6.65(3) 8.43(4) 13.85(6) 

U E /Ne -1.511(1) -1.695(1) -1.778(2) -1.898(1) -2.116(1) -2.396 -2.754 -3.214 -3.817 

T* =1.0 



Z 7.58(2) 8.96(2) 6.81(3) 6.88(2) 7.20(3) 7.85(3) 9.09(3) 11.74(4) 18.23(4) 

U E /Ne -1.465(1) -1.647(1) -1.673(1) -1.849(1) -2.078 -2.368 -2.734 -3.204 -3.813 

T* =1.5 



Z 8.63(2) 10.25(2) 7.96(2) 8.33(2) 9.15(2) 10.42(2) 12.45(3) 16.23(3) 24.36(4) 

U E /Ne -1.402(1) -1.587(1) -1.603(1) -1.787 -2.029 -2.331 -2.710 -3.191 -3.809 

T* =2.0 



Z 9.22(2) 10.82(2) 8.46(2) 9.11(2) 10.20(2) 11.71(2) 14.16(2) 18.51(3) 27.37(3) 

U E /Ne -1.376(1) -1.558 -1.565(1) -1.759 -2.006 -2.313 -2.698 -3.186 -3.807 

T* =3.0 



Z 9.73(2) 11.58(2) 9.03(2) 9.86(2) 11.20(2) 13.09(2) 15.89(2) 20.77(3) 30.41(3) 

U E /Ne -1.347(1) -1.536 -1.533 -1.732 -1.984 -2.297 -2.687 -3.180 -3.805 
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TABLE IV: As in Table Q] for 7 = 36. 



p* 0.90 0.95 1.00 1.05 1.10 1.15 1.20 1.25 1.30 
T* =0.6 

Z 6.17(2) 7.20(2) 7.24(6) 5.95(3) 5.86(3) 5.79(3) 5.47(3) 4.88(4) 4.11(6) 

U E /Ne -1.061(2) -1.194(1) -1.284(5) -1.319(2) -1.460(2) -1.648(2) -1.890(1) -2.233(1) -2.753 

T* =0.8 



Z 7.58(2) 8.89(3) 7.72(4) 7.61(3) 8.00(3) 8.47(3) 9.03(3) 9.91(5) 12.05(6) 

U E /Ne -0.960(1) -1.097(1) -1.100(2) -1.192(2) -1.341(1) -1.536(1) -1.794(1) -2.160 -2.716 

T* =1.0 



Z 8.42(2) 9.80(2) 8.11(2) 8.53(2) 9.07(3) 9.95(3) 11.07(3) 12.92(4) 16.78(6) 

U E /Ne -0.907(1) -1.033(1) -1.012(1) -1.126(1) -1.272(1) -1.471(1) -1.739(1) -2.119 -2.695 

T* =1.5 



Z 9.32(2) 11.02(3) 8.81(2) 9.49(2) 10.54(2) 11.90(2) 13.89(3) 16.97(4) 23.19(4) 

U E /Ne -0.838(1) -0.970(1) -0.919(1) -1.034(1) -1.188 -1.393 -1.673 -2.070 -2.669 

T* =2.0 



Z 9.65(2) 11.45(2) 9.22(2) 10.01(2) 11.21(2) 12.85(3) 15.28(2) 19.09(3) 26.46(4) 

U E /Ne -0.802 -0.930(1) -0.877(1) -0.992 -1.149 -1.357 -1.642 -2.047 -2.657 

T* =3.0 



Z 10.03(2) 11.59(6) 9.56(2) 10.51(2) 11.89(2) 13.82(2) 16.66(2) 21.19(3) 29.87(4) 

U E /Ne -0.771 -0.883(3) -0.838 -0.954 -1.112 -1.321 -1.611 -2.026 -2.647 
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List of Figures 



FIG. 1: Comparison between MC simulations (filled circles) and perturbation theory for the excess 
energy U E /Ne of Sutherland solids for T* = 0.6, 0.8, 1.0, 1.5, 2.0, and 3.0 from the bottom upwards. 
Open circles are the results from the MC perturbation theory as obtained from eq. (I20p . MC and 
MC-P results are indistinguishable at the scale of the figure in most cases. Dotted curves are the 
results from the first-order Barker-Henderson perturbation theory. Dashed and continuous curves 
are the results from the second-order Barker-Henderson perturbation theory using the macroscopic 
and local compressibility approximations, respectively. In most cases, the three curves are nearly 
indistinguishable to each other at the scale of the figure. For each temperature the curves, and the 
corresponding sets of data, have been displaced upwards for clarity by 1 unit with respect to those 
corresponding to the temperature immediately below. 

FIG. 2: Comparison between MC simulations (filled circles) and perturbation theory for the com- 
pressibility factor Z = pV/NksT of Sutherland solids for T* = 0.6, 0.8, 1.0, 1.5, 2.0, and 3.0 from 
the bottom upwards. Open circles are the results from the MC perturbation theory as obtained 
from eq. (|2ip . Again, MC and MC-P results are nearly indistinguishable to each other. The curves 
have the same meaning as in fig. [TJ For each temperature the curves, and the corresponding sets 
of data, have been displaced upwards for clarity by 5 units with respect to those corresponding to 
the temperature immediately below. 
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